Hello world
plot(runif(100))

dat<-read.csv("_dat/richness-data.csv")
head(dat)
## X ele_meters richness ln_ele ln_richness
## 1 1 1383 20.51757 7.232010 3.021282
## 2 2 764 15.61994 6.638568 2.748548
## 3 3 467 14.53794 6.146329 2.676762
## 4 4 457 10.15791 6.124683 2.318253
## 5 5 1779 17.97125 7.483807 2.888773
## 6 6 1415 13.06136 7.254885 2.569658
plot(richness~ele_meters,dat)

fit<-lm(ln_richness~ln_ele,dat)
fit
##
## Call:
## lm(formula = ln_richness ~ ln_ele, data = dat)
##
## Coefficients:
## (Intercept) ln_ele
## -0.02814 0.40713